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Abstract 

Let S'^ denote the unit sphere in the Euclidean space {d > 1). Let A be a 

natural number {N > 2), and let ujn ■= • • •, xat} be a collection of N distinct 

points on on which the minimal Riesz s-energy is attained. In this paper, we show 
that the points xi,... ,X]\f are well-separated for the cases d — 1 < s < d. 


1 Introduction. 


Let 5'^ denote the unit sphere in the Euclidean space {d> 1). Let A be a natural 
number (A > 2), and let ujn := {xi,... ,xn} be a collection of A distinct points on 
The Riesz s-energy (s > 0) associated with lon, Es{^n), is defined by 


Es{ujn) ■= 


Xj X -j 


E 

I—y—|, 


if s > 0, 
if s = 0. 


Here | • | denotes the Euclidean norm. We use Ss{S'^,N) to denote the A-point 
minimal s-energy over S'^ defined by 


Ts(5‘^,A):= min Es{ujn), (1-1) 

ujnCS<^ 
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where the infimum is taken over all A/^-point subsets of S'^. If lon C S'^ is such that 

Es{u;n)=£s{S^N), 

then Wat is called a minimal s-energy configuration, and the points in cuat are called 
minimal s-energy points, or simply minimal energy points if the linkage to the pa¬ 
rameter s is well-understood in a certain context. It is obvious that minimal s-energy 
configurations exist. Also, if ujn is a minimal s-energy configuration, and if p is a 
metric space isometry from to 5^, then the image of ujn under p, p{uj]\f), is also a 
minimal s-energy configuration. Minimal (d — l)-energy points are often referred to as 
“Fekete points”; see [3]. The determination of minimal s-energy configurations and the 
corresponding minimal s-energy on S'^ and other manifolds is an important problem 
that has applications in many subjects including physics, chemistry, computer science, 
and mathematics. For further background regarding this problem and its applications, 
we refer readers to the expository papers by Hardin and Saff na, and by Saff and 
Kuijlaars m- The papers M, i, m, m, and m and the references therein also 
contain valuable pertinent information. 

The determination of the distribution of minimal s-energy points in the cases d > 2 
turns out to be rather elusive. It is, however, generally expected that these points are 
“well-separated” in the sense that there exists a positive constant depending only 
on d and s, such that 

m.m\xi - Xj\ > (1.2) 

Dahlberg |3j proved that Fekete points are well-separated. (Dahlberg actually estab¬ 
lished the well-separatedness of Fekete points on every compact C^’^-surface in 
Kuijlaars and Saff m proved that minimal s-energy points are well-separated for the 
cases s > d. There have been a series of more quantitative results regarding the case 
d = 2, s = 0, which corresponds to the minimal logarithmic energy on Rakhmanov, 
Saff, and Zhou m first showed that 3/5 is a lower bound for the constant A 2 ^o in in¬ 
equality dca). Dubickas [SI refined Rakhmanov, Saff, and Zhou’s method, and showed 
that 7/4 is a lower bound for the constant ^ 2,0 in inequality (11.21) . Using potential 
theory and stereographical projection techniques, Dragnev 0 established an appealing 
lower bound for the minimum separation of the minimal logarithmic energy points on 
52 to be 2(iV- 1)-V2. 

In this paper, we show that the minimal s-energy points are well-separated for 
the cases d — 1 < s < d. Note that the case s = d — 1 is already covered by the 
aforementioned Dahlberg’s result. While the parameter s is restricted in the range 
d — 1 < s < d, the outcome s = 0 can only occur when d = 1, putting the problem 
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on the unit circle S^. It is shown by Gotz m that the Nth roots of unity and their 
rotations are the only minimal s-energy configurations on ; see also m- Therefore in 
this paper, we can use the tacit assumption that s > 0. Our proof entails comparing the 
Riesz s-potentials on the slightly larger sphere of radius 1 + of two probability 

measures: the rotationally invariant probability measure on and the normalized 

counting measure on a minimal s-energy configuration a; at. The Riesz s-potential of 
the rotationally invariant probability measure on S'^ can be expressed in closed form in 
terms of the Gauss hypergeometric functions 2-^1 (oj b; c; z). The crux of our argument 
is an application of the principle of domination for a-superharmonic functions; see 
m This paper is organized as follows. In Section 2, we introduce the necessary 
notations and terminologies. Also in Section 2, we list a few formulas pertaining to 
the hypergeometric functions that we will use in our proofs. In Section 3, we state and 
prove our main result. 


2 Notation and terminology. 


Given a minimal s-energy configuration lon, we use vw to denote the normalized count¬ 
ing measure on a; at, i.e., 

N 

VN ■= 

i=i 

where 6xj denotes the unit point mass at Xj. The rotationally invariant probability 
measure on S"^ is denoted by 11 . For a cj-finite positive Borel measure A supported on 
a compact subset K of we define its Riesz s-potential (s > 0) by 


Us{x) ■= [ \x-y\ *dA(y). 

JK 


Note that Ug may take the extended value 00 on some subsets of In this paper, 

K is either or wat- In the next section, we will show that the Riesz s-potential 
of y, can be expressed in closed form in terms of the Gauss hypergeometric functions 
2 F 1 (a, b] c; z) defined by 


00 

2 Fi(a,b;c-,z) := ^ 
n=0 


(a)n(b)n 

(c)n n!’ 


where (a)„ is the Pochhammer symbol defined by 


{a)n ■= 


if n = 0, 


1 , 

a(a -|- 1) • • • (a -|- n — 1) 


if n > 1. 
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There are several sources that provide essential properties of the hyper geometric func¬ 
tions 2 ^ 1 ( 0 , b; c] z) such as Abramowitz et al P, and Andrews et al [2]. We will be 
using a few basic formulas pertaining to the hyper geometric functions 2-^1 (o, c; z), 
which can all be found in We quote them here for easy reference. Under the 
conditions Re (c) > Re (6) > 0, the following formula holds true: 

[\l - zu)-'^u’’-\l - uy-’^-^du = 2 Fi{a, b- c; z). (2.1) 

Jq r(c) 

The above integral represents an analytic function in the z-plane cut along the real 
axis from 1 to 00 . Formula (EH) is often called Euler’s integral representation for the 
hypergeometric function 2 F 1 . If Re(c — a — b) < 0, then 


,. 2 ^ 1 ( 0 , b; c; z) _ r(c)r(a + b - c) 

. (1 - zy-<^->^ ~ r(a)r(6) 

If Re(c — a — b) > 0, then 


( 2 . 2 ) 


2 Fi(a,b; c; 1) 


r(c)r(c — a — b) 
r(c — a)r(c — b) 


(2.3) 


The following derivative formula can be easily proved by term-by-term differentiation 


in a suitable domain: 


d_ 

dz 


2 Fi{a,b]c; z) 


ah 

c 


2 Ti(a-|- l,b l;c-|- l;^:). 


(2.4) 


3 Main Result and Proofs. 

On various occasions, we use Cd,s to denote some unspecified positive constants, de¬ 
pending only on d and s. The exact values of may be different from proof to 
proof. In the same proof, however, for clarity we use different notations for different 
constants, namely C'^ and if necessary. Although in the current paper, we do not 
strive to estimate these constants, closed forms for them can be obtained with some 
devoted calculations. We will be using the notations former 

denotes a “double sum”, excluding only those terms with i = j. The latter denotes a 
single sum in which i is fixed, and the summation is done only on j. 

The following lemma is well-known; see e.g. 1111 - 

Lemma 1. For d — 1 < s < d, there exists a positive constant independent of N, 
sueh that for any N distinct points xi,... ,X]\r on we have 

Y, \xi - > id,sN^ - 


(3.1) 






where 
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^d,s 



y\ ^dn{x)dfj,{y) 


r{{d + i)/2)r{d-s) 

T{{d- s + l)/2)T{d- s/2)' 


By inspecting the pertinent proof in m, we find that a quantitative estimate of 
the constant Cd,s in the above lemma is possible. The determination of asymptotically 
sharp values for these constants has received much attention in the literature; see |3|, 
m, m, and references therein. 

Lemma 2. ForO < s < d, there exists a positive constant Cd,s independent of N, such 
that 


Uf^{x) > ^d,s - |x| = 1. 


(3.2) 


Proof: Since vn is the normalized counting measure of a minimal s-energy 

configuration lon = {xi,..., xat}, we have for each fixed i, 1 < i < N, 


Y, \x-xj\ * > ^ 


\ •} 




Summing over i and using Lemma ^ we get 


N 


{N — 1) \x — X 




(3.3) 


Dividing by N{N — 1) we get the desired estimate. 


Lemma 3. For each fixed s > 0, the potential Us of the measure y is a radial function, 
and has the explicit expression in terms of the hypergeometric function: 


' s d 47? 

U^{x) = (7? + l)-2Ti(-,-;d; 


|x| = 7? / 1. 


Proof: For x G 


2’ 2’“’ (7? + 1)V ’ 

, |x| / 1, we have 

U^{x)= f \x-y\~^dy,{y). 

Jsd 

Let x,y G with \x\ = 7? and \y\ = 1. Denote the angle between the two vectors x 
and y by 6. Then cos0 = (-§,2/). By the law of cosine, \x — y\‘^ = Ef + 1 — 2R{^,y). 
Thus by using the Funk-Hecke formula; see Muller [m, we have 


ISd 


\x-y\ ®d/i(y) = [ (7?^ + 1 - 27?(^,2/)) ^^^dp{y) 

Jsd- E 

J (7?2 + l-27?7)-'*/2(l_i2)(d-2)/2^^^ 


Isd 

^d-1 


^d J-l 
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where I'd denotes the surface area of S'^. Using the substitution 2u = t + l and Euler’s 
integral representation of the hyper geometric function 2 -E 1 , we have 


U^{x) 




T^d-l 

i^d 

^d-i 

Vd 



AR 


{R + lf 


s d , 

2 ’ 2 ’ ’ 


-s/2 




4R \ 

(imyj- 




To simplify, we use the formula (see ini) 


= 


27r('^+i)/2 

r(^) ’ 


and then the formula (see pQ) 


r(2z) = (27r)-^/^2^^-^/^r(z)r(z + 1 / 2 ), 


(3.4) 


with z = d/2. We have 


U,^(x) 


2'^-\R-hl)-^ 
(i? + l)-^Fi 


r(^) r^(l) 
m' ' 

(sd 4R ' 
V2’2’"^’ {R + l)\ 


2 ’ 2 


;d; 


m \ 

w+w) 


The following two special cases of Lemma |21 are worth noting. Firstly, when 0 < 
s < d, the potential Us{x) is well defined for |x| = 1. In fact, a simple application of 
the Lebesgue Dominated Convergence Theorem yields 


lim U^{x) 


lim / \x-y\ ^^d/iiy) ='yd,s- 


On the other hand, when 0 < s < d, the hyper geometric function 2 -E 1 (|,|;d; 2 ;) is 
continuous at 2 ; = 1. Using equation (lOl) . we have 


lim {R+ 1 )-^ 2 Fi 
R^ l 


(sd AR \ 
V2’2’“^’ {R+lfJ 


2 *2Ti(|,^;d;l) 

r(d)r(^) 

r(d/ 2 )r(d-s/ 2 )' 


Thus, we have 


'JdjS — 2 


r(d)r(^) 
r(d/2)r(d- s/2)' 


(3.5) 


Due to the importance of the constant 'yd,s in this paper, we feel reassured that 
we are able to verify equation (13.511 directly, and we share the reassurance with the 
readers. Using equation (|S3I) with 2 ; = d/2, we write 


r(d) = (27r)-^/22rf-i/2p(^/2)r((d + l)/ 2 ) 
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we have 


r(d)r(^) ^ (27r)-v^2"-v^r(d/2)r((d + i)/ 2 )r(^) 

T{d/2)T{d- s/2) T{d/2)T{d- s/2) 

(27r)-i/22<i-s-i/2r(^^)r('^+i)r((d + l)/2) 

r{^=f^)T{d- s/2) 

r(d-g)r((d + i)/ 2 ) 

r(iLi|±i)r(d-s/2) ■ 

Here in the last step, we have used equation (|S3I) again with z = {d — s)/2. Secondly, 
when d = 2, the potential Us{x) has the elementary form: 


U^{x) 


1 (l + i?)2-^ - |i?- 1|2-^ 

m 2 - s 


x\ = R, 


s ^ 2. 


Lemma 4. Assume d — 1 < s < d. Then there exists a positive constant Cd^sj inde¬ 
pendent of N, such that 


Us^{x) > 

for all X G with |x| = 1 + 


Proof: We first note that with R]sf,d := 1 + iV we have 


{RM,d + 1)“" = 2-*(l - + o{N-^/'^). (3.6) 

We now estimate 2-^1 j f! ^^5 (^^ 4 ^)' using the Fundamental Theorem of Cal¬ 
culus and then equation (ITU) . we have 


2-ci 


•s d , 4i?jv,d 
2 ’ 2 ’"'’ {RN,d + ir 


= 2 F 1 


s d 
2’2 


;d; 1 - 


2-ci 


s d , , ^ 

-,-;d;l) - 2 T 1 


r(d)r( 


d—s ^ 
2 > 


r(d/2)r(d - s/2) 4 74fljv,d/(Rjv.d+i)2 


f 

JAR. 


2 FI 


sd 4i?Ar,d 
2’2’^’ {RN,d + ^? 


— 1, — + 1; d + 1; z ) dz. 


We use equation (EH to estimate the above integral. For any given e > 0, we have, 
for N sufficiently large, that 


2-^1 (f + 1) f + 1; d + 1; 2 ;) 

(1 _ _j)(d-s-2)/2 



< e. 


Z G [4i?Ar,d/(i?Ar,d + 1)^, 1], 


F(d+l)F(((i-s)/2 + l) 
F(d/2 + l)F(s/2 + 1) 


where 
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This implies 


L 


2-^1 


4-RiV,d/{-RjV,d + l)^ 


- + I, - + 1] d + I] z ] dz 


< iPd,s + e) C (1 - z)^^-^-^y^dz 

j 4-Rjv,d/(i?JV,d + l)^ 


— {Pd,s + c)‘ 


Rna — 1 


d — S \ RN,d + 1 

OS —d+1 


d—s 


(3.7) 


Combining Lemma 3 and the above two estimates (ESI) and EIZI), and noting that 
1/d > 1 — s/d, we have for x with |x| = RN,d, 


U^ix) 


2-"(l 


£iv-Vrf) 


T{d/2)T{d- s/2) 



+ o{N-^+^/^), 


which gives the desired result of Lemma 0] ■ 

In the discussion that follows, we will need the notion of “a-superharmonic func¬ 
tions” and the principal of domination of a-superharmonic functions. These topics 
can be found in Landkof ESI- The definition of a-superharmonic functions is tech¬ 
nical. Upon checking the pertinent material in Landkof ESI, one finds the relation 
d + 1 — a = s between the parameter a used in Landkof ESI and the parameter s we 
use here. Thus the requirement d — 1 < s < d translates into 1 < a < 2 in Landkof ESI- 
Furthermore, what is meant in Landkof m by an a-superharmonic function is in fact 
a (d -|- 1 — s)-superharmonic function here in our context. We choose not to use the 
phrase (d-|-1 — s)-superharmonic function because we feel that notion a-superharmonic 
function has been coined in the mathematical literature. In Chapter 1, Section 5 of 
ESI, it is proved that for a u-finite Borel measure A supported on a compact subset 
of its potential Ug (0 < s < d) is an a-superharmonic function. Furthermore, 

Theorem 1.29 in ESI, aside from some changes in notation, states the following result: 


Theorem 5. Suppose X is a a-finite positive Borel measure supported on a compact 
subset ofW^^^ whose potential is finite X-almost everywhere, and that f{x) is an 
a-superharmonic function. If the inequality 


U^x) < fix) 


holds X-almost everywhere, then it holds everywhere in 

Theorem 13 is often called the “principal of domination” for a-superharmonic func¬ 


tions. 
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Lemma 6. Assume d — 1 < s < d. Then there exists a constant Cd,s, independent of 
N, such that for all x G with |x| = 1 + 

C/r (x) > ld,s - Cd,sN-^+^/^. (3.8) 

Proof: By Lemma [21 we have 

Urix) > ^d,s - Cd,sN-^+^/<^, |x| = 1. 

Note that Us{x) = '^d,si x G so we can rewrite the above inequality as 

Urix) > U^{x){l - |x| = 1. 

Both measures /x and vn are supported on S'^. Since is an a-superharmonic 
function, we can therefore use Theorem 0 to obtain for N sufficiently large, 

Uf^{x) > U^{x){l - x G 

We then use Lemma 0] to get for |x| = 1 + 

Uf^ix) > i^d,s-ClN-^+^/^){l-C'd,,N-^+^/<^), 

which yields the desired estimate. ■ 

Lemma 7. For every fixed i,{l < i < N), we have, for 0 < s < d, 

^ \xi - Xj\~^ < 'yd,s- 

Proof: For every fixed i, {1 < i < N), since the function x 'Ylj-jjki I® “ Xj\~^ 

reaches its minimum at Xj on S"^, we have 

\xi — Xj\~^ < N~^ |x — Xj\~^, x G S’^. 

Integrating both sides of the above inequality on against the probability measure 
/i(x) yields: 

X] I®* “ /rf X] I® “ Xj\~^d^{x) = ^^ ^ id^s, 

and therefore the desired inequality follows. ■ 

Theorem 8. Assume d — 1 < s < d. The minimal s-energy points are well-separated, 
i.e., there exists a constant Ad^s > 0, independent of N, such that 

min \xi — Xj\ > Ad^s^~~^^'^■ 
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^d,s ^ 1^*0 ®io I 


Proof: Let ujn be a minimal s-energy configuration, and let Xi^ , Xj^ be two points 
in Wat such that 

l.r,-. — .r,'. I =_ 

i¥=j 

Using Lemma Q we have 

> iv-> 

j-j¥=io 

= ^ \xif^-Xj\~\ (3.9) 

j-j¥=io,jo 

Take x := (1 + N~^/^)xiQ. Then |xjg — Xj\ < |x — Xj\ for every j, and so by (EH and 
Lemma ini we have 

7-1 


I I ^ I ^io ^ jo I 


1(1,3-N ^\xiQ-Xj^\ ^ > N 


E I- 

iuV*ojo 

= {x) - N-^\x - x^ 

> ld,3-Cd,sN-^+^/‘^ -N-^\x-. 
Since |x — Xjpl = and |x — x^qI > it follows that 


*01 ^ I ® ®io I 


'iol ^ ^ V - a^iol *■ 


ld,s ^ l^io 


> ld,s-{Cd,s + 2)N-^+^/<^, 

which implies that for some constant 


XlQ X 


30 


> Ad,sN-^/<^. 


If one follows the trail of the constants throughout the proofs, then one is able to 
quantitatively estimate the value of the constant in Theorem |H1 However, such a 
process leads to the piling-up of many Gamma function values, among others. There 
seems to be no obvious way to simplify them. Just for curiosity, we numerically esti¬ 
mated the constant Ad,s in Theorem |H1 for the case d = 2, and s = 1. Our numerical 
result yields A 2 J > 0.8709, putting our estimate of the minimum separation of the 
corresponding minimal energy points on at 0.8709/\/A/'. The result of Habicht 
and Van der Waerden EH for best packing asserts that the maximum diameter of N 
non-overlapping congruent circles on is asymptotically 


/ 87r \ 1 

[TsJ 7W 


3.809 


1 


^/n' 


For more information on best packing on S^, we also refer readers to n, m, and the 
references therein. 
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